Astronomical Journal, in press (April 2001) 

Preprint typeset using LAT^X style emulateapj v. 20/04/00 



O 
O 

o 

O 

Q 

oo 



oo 
m 

(N 

o 

O 

Or 

6 

in 
c3 



DYNAMICAL FAMILY PROPERTIES AND DARK HALO SCALING RELATIONS OF GIANT 

ELLIPTICAL GALAXIES 

Ortwin Gerhard, Andi Kronawitter 

Astronomisches Institut, Universitat Basel, Venusstrasse 7, CH-4102 Binningen, Switzerland; 
gerhard@astro.unibas.ch, krona@unibas.ch 

R.P. Saglia, Ralf Bender 

Institut fur Astronomie und Astrophysik, ScheinerstraBe 1, D-81679 Munich, Germany; 
saglia@usm.uni-muenchen.de, bender@usm.uni-muenchen.de 
Astronomical Journal, in press (April 2001) 

ABSTRACT 

Based on a uniform dynamical analysis of the line-profile shapes of 21 mostly luminous, slowly rotating, 
and nearly round elliptical galaxies, we have investigated the dynamical family relations and dark halo 
properties of ellipticals. Our results include: (i) The circular velocity curves (CVCs) of elliptical galaxies 
are flat to within ~ 10% for R > 0.2R e . (ii) Most ellipticals are moderately radially anisotropic; their 
dynamical structure is surprisingly uniform, (iii) Elliptical galaxies follow a Tully-Fisher (TF) relation 
with marginally shallower slope than spiral galaxies, and -y™ ax ~ 300 km s -1 for an L* B galaxy. At given 
circular velocity, they are ~ 1 mag fainter in B and ~ 0.6 mag in R, and appear to have slightly lower 
baryonic mass than spirals, even for the maximum M/ Lb allowed by the kinematics, (iv) The luminosity 
dependence of M/Lb indicated by the tilt of the Fundamental Plane (FP) is confirmed. The tilt of the 
FP is not caused by dynamical or photometric non-homology, although the latter might influence the 
slope of M/L versus L. It can also not be due only to an increasing dark matter fraction with L for 
the range of IMF currently discussed. It is, however, consistent with stellar population models based 
on published metallicities and ages. The main driver is therefore probably metallicity, and a secondary 
population effect is needed to explain the K-band tilt, (v) These results make it likely that elliptical 
galaxies have nearly maximal M/Lb (minimal halos). (vi) Despite the uniformly flat CVCs, there is 
a spread in the luminous to dark matter ratio and in cumulative M/Ls(r). Some galaxies have no 
indication for dark matter within 2R e , whereas for others we obtain local M/LbS of 20-30 at 2R e . (vii) 
In models with maximum stellar mass, the dark matter contributes ~ 10 — 40% of the mass within R e . 
Equal interior mass of dark and luminous matter is predicted at ~ 2 — 4R e . (viii) Even in these maximum 
stellar mass models, the halo core densities and phase-space densities are at least ~ 25 times larger and 
the halo core radii ~ 4 times smaller than in spiral galaxies of the same circular velocity. The increase 
in M/L sets in at ~ 10 times larger acceleration than in spirals. This could imply that elliptical galaxy 
halos collapsed at high redshift or that some of the dark matter in ellipticals might be baryonic. 

Subject headings: galaxies: elliptical and lenticular - galaxies: kinematics and dynamics - galaxies: 
stellar content - galaxies: halos - galaxies: formation - cosmology: dark matter 



1. INTRODUCTION 

Hierarchical theories of galaxy formation predict that el- 
liptical galaxies should be surrounded by dark matter ha- 
los. The study of these halos has been difficult, however, 
because of the lack of suitable and easily interpreted trac- 
ers such as the HI rotation curves in spiral galaxies. Re- 
cent work in several fields, however, leaves no doubt about 
the existence of dark matter in ellipticals: X-ray data on 
their hot gas athmospheres implies that dark halos in el- 
lipticals are ubiquitous and that the mass-to-light ratios 
are M/L ~ 100 on scales of ~ 100 kpc (Mushotzky et al. 
(1994), Matsushita et al. (1998), Loewenstein & White 
(1999)). Gravitational lensing studies (Kochanek (1995), 
Keeton, Kochanek & Falco (1998), Griffiths et al. (1998)) 
show evidence for large M/L in lens elliptical galaxies, and 
stellar-dynamical studies based on absorption line profile 
shapes have given strong constraints on the mass distribu- 
tions to ~ 2R e (Rix et al. (1997), Gerhard et al. (1998), 
Saglia et al. (2000a)), with a small to moderate dark mat- 
ter fraction inside 2R e . 

Despite of this progress, the detailed mass distributions 



in elliptical galaxies and their variations with luminosity 
remain largely unknown. Has the luminous matter segre- 
gated dissipatively in the halo potential? Is there a "con- 
spiracy" between luminous and dark matter to produce 
a flat rotation curve, like in spiral galaxies? How do the 
mass-to-light ratio, the slope of the circular velocity curve, 
or the orbital anisotropy scale with luminosity? Is the tilt 
of the Fundamental Plane simply related to a variation of 
the dynamical M/L with LI Do elliptical galaxies follow 
a Tully-Fisher relation? How do the scale radii and den- 
sities of elliptical galaxy halos compare to those of spiral 
galaxies? 

The purpose of this paper is to address some of these 
questions on the basis of a new dynamical study by Kro- 
nawitter et al. ((2000), hereafter K+2000), who analyzed 
the dynamical structure and mass distribution for a sample 
of 21 bright elliptical galaxies. Continuing from previous 
work by our group (Gerhard et al. (1998), Saglia et al. 
(2000a)), these authors modelled the line profile shapes 
of in total 17 E0/E1 and 4 E2 galaxies for which kine- 
matic data including line profile information were avail- 
able from their own observations or from the literature. 
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The dynamical structure of these galaxies turned out to 
be remarkably uniform. Most galaxies require moderate 
radial anisotropy in their main bodies (at ~ 0.5R e ). Their 
circular velocity curves are all consistent with being flat 
outside ~ 0.2R e . The M/L ratio profiles begin to rise at 
around 0.5 — 2R e and are consistent with X-ray and other 
data where available, although from the kinematic data 
alone constant M/L models can only be ruled out at 95% 
confidence in a few galaxies. 

This sample provides a new and much improved basis 
for investigating the dynamical family properties of ellip- 
tical galaxies, which is the subject of the present study. In 
Section 2, we analyze the unexpectedly uniform dynami- 
cal structure of these elliptical galaxies. In Section 3, we 
investigate the dependence on luminosity, discussing the 
Faber-Jackson, Tully-Fishcr and Fundamental Plane rela- 
tions. In Section 4, we relate the dynamical mass-to-light 
ratios to the stellar population properties. In Section 5 we 
discuss the structure of the dark halos of these ellipticals. 
Our conclusions are summarized in Section 6. 

2. DYNAMICAL STRUCTURE 

The elliptical galaxies analyzed by Kronawitter et al. 
((2000), K+2000) divide in two subsamples, one with new 
extended kinematic data, reaching typically to ~ 2R e 
('EK sample', the data are from Kronawitter et al. and 
from several other sources referenced there) , and one with 
older and less extended kinematic measurements ('BSG 
sample'; this is a subsample from Bender, Saglia & Ger- 
hard (1994)). Based on these data and mostly published 
photometry, K+2000 constructed non-parametric spheri- 
cal models from which circular velocity curves, radial pro- 
files of mass-to-light ratio, and anisotropy profiles for these 
galaxies were derived, including confidence ranges. 

The galaxies were selected to rotate slowly if at all and 
to be as round as possible on the sky. They are lumi- 
nous elliptical galaxies (Mb — —21 ± 2) 1 . The expected 
mean intrinsic short-to-long axis ratio for such a sample of 
luminous ellipticals is < c/a>= 0.79. The mean system- 
atic effects arising from the use of spherical models and 
the possible presence of small embedded face-on disks are 
small for the sample as a whole, but may be non-negligible 
in individual cases (see K+2000, §5.1). 

2.1. Circular velocity curves 

Circular velocity curves (CVCs) for all galaxies in the 
sample are shown in Figure 1, in three bins roughly or- 
dered by luminosity. CVCs normalized by the respective 
maximum circular velocity are shown in Figure 2 sepa- 
rately for the two subsamples. The plotted curves corre- 
spond to the "best" models of K+2000, which are taken 
from the central region of their 95% confidence interval for 
each galaxy, respectively. Based on dynamical models near 
the boundaries of the confidence interval, the typical un- 
certainty in the outermost circular velocity is ±(10-15)%. 
The expected mean systematic error from flattening along 
the line-of-sight is smaller; cf. §5.1 of K+2000. 

The most striking result from these diagrams is that at 
the ~ 10% level all CVCs are flat outside R/R e ~ 0.2. 
This result is most significant for the galaxies with the ex- 
tended data, while for many galaxies from the BSG sample 

throughout this paper we use a Hubble constant Ho = 65kms~ 



the radial extent of the data is insufficient to show clear 
trends. However, in cases where X-ray data are available 
(NGC 4472, 4486, 4636) the mass profiles of the "best" 
models approximately match those from the X-ray analy- 
sis even for those galaxies (see K+2000). 

This result is illustrated further by Figure 3, which 
shows the derived ratio v c (R max ) / v™ aK for all galaxies of 
the EK-sample. Here v c (R max ) is the circular velocity at 
the radius of the last kinematic data point, and ^™ ax is 
the maximum circular velocity in the respective "best" 
model. For NGC 315 v c (0.6R e ) was used instead of < ax . 
The error-bars plotted correspond to the 95% confidence 
range for w c (ii max ), compared to which the uncertainty 
in v™ ax can be neglected. Most of these galaxies have 
«c(^max)/f c max ^ 0.9 - 1.0 with a median at 0.94, and 
95% confidence ranges ~ ±0.1. 

The galaxies in the EK sample appear to show a "bi- 
modal" distribution of CVC shapes in Figs. 1, 2: For one 
group of galaxies the "best" model CVC has a peak near 
0.3i? e and then falls slightly until it becomes flat at ~ lR e 
(NGC 1399, NGC 3379, NGC 5846 and NGC 6703, the 
latter has the largest drop). In the other group the CVC 
rises rapidly until ~ 0.2R e and reaches a peak at ~ lR e 
after which it remains nearly flat (NGC 2434, NGC 7145, 
NGC 7192, NGC 7507). NGC 4374 and NGC 7626 may 
be cases from the first group where the flat part is not 
yet seen in the data. The rise of the CVC of NGC 315 
seen near i? max in Fig. 1 may not be real; a flat CVC 
appears also consistent with the kinematic data (see dis- 
cussion in K+2000). The difference between these CVC 
shapes is about a ~ 2a result when comparing the model 
confidence bands for EK galaxies of either type, but this 
estimate does not include possible systematic effects from 
the use of spherical models. If the differences in CVC 
shapes are real, they could reflect small variations in the 
degree of dissipation and mass segregation of the baryonic 
component during the formation process. In any case no 
clear trend with galaxy luminosity is seen. 



2.2. Anisotropy 

Figure 4 shows the radial profiles of the anisotropy pa- 
rameter (3 = 1 — of /of for all galaxies. Here oy and a t 
are the intrinsic radial and one-dimensional tangential ve- 
locity dispersions, and (3 = 1,0, — oo for completely radial, 
isotropic, and circular orbit distributions, respectively. Be- 
cause it is a deprojected quantity, the anisotropy derived 
from our models is considerably more uncertain than the 
total mass M(r); see Gerhard et al. ((1998)) for further 
discussion. The uncertainties are particularly large near 
the outer boundary of the data, where the range of [3 val- 
ues in models corresponding to the allowed potentials is 
usually 0.2-0.5. For the inner profiles, the allowed range 
is typically 0.1-0.2, sometimes 0.3. Nonetheless, Figure 
4 shows a clear trend in that almost all galaxies are ra- 
dially anisotropic in the inner (best-constrained) regions, 
with the peak of the anisotropy often near 0.2 R e with 
values of (3 = 0.2. ..0.4. There are two cases in the EK 
sample which are exceptional in that they are consistent 
with isotropy over the whole radial range. These two are 
NGC 315 (thick solid line), a cD galaxy, and NGC 7626 

unless explicitly noted otherwise. 
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Fig. 1. — The "best model" circular velocity curves of all galaxies from the K+2000 sample plotted as a function of radius scaled 
by the effective radius R e . The panels are roughly ordered by luminosity. 
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Fig. 2. — Same circular velocity curves, normalized by the maximum circular velocity. The upper panel now shows the galaxies 
from the EK subsample of K+2000, the lower panel those from the BSG subsample. The extended curve in the lower panel is for 
the compact elliptical NGC 4486B. 



(thick long dashed line), which is classified as E pec. De- 
spite their less extended and lower S/N data, the galaxies 
from the BSG sample show a similar trend, again with 
two exceptions: NGC 4636, which we do not consider a 



very reliable case, and NGC 4486B, which is exceptional 
in almost all respects. 

To get a more robust estimate of the anisotropy we 
have averaged the anisotropy profile between 0.1 R e and 
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Fig. 3. — The ratio w c (-Rmax)/w™ ax for all galaxies from the EK-sample, showing the gradient in the circular velocity curves. The 
error-bars plotted correspond to the range of v c (i? ma x) from the 95%-confidence models of K+2000. 



min(i? e , -Rmax), where J? ma x is the radius of the outermost 
kinematic data point. These averaged anisotropies are 
plotted in Figure 5 versus the maximum circular veloc- 
ity from Fig. 1. In a few cases this is the circular velocity 
at the outer boundary of the modelled range. Most galax- 
ies with extended data show a clear maximum of the CVC 
whereas the BSG galaxies have CVCs which may still rise 
outwards. For these the true maximum probably lies out- 
side the range of the kinematic data. The error-bars plot- 
ted for (3 show the standard deviation from the mean over 
the radial range used in the averaging. The errors plotted 
for v c are equal to ± one half the separation of the extreme 
models in the confidence interval, near the radius where 
the best model has its maximum. The figure shows again 
that all galaxies but two (NGC 4486B and NGC 4636) 
have a mean radial anisotropy in the range from (3 ~ to 
(3 ~ 0.35, and also that there is no dependence on circular 
velocity (or luminosity). 

Radial anisotropy has also been inferred from three- 
integral axisymmetric models for several flattened ellip- 
tical galaxies: NGC 1600 (E3.5, Matthias & Gerhard 
(1999)), NGC 2300 (E2, Kaeppeli (1999)), NGC 2320 
(E3.5, Cretton, Rix & de Zeeuw (2000), NGC 3379 (El, 
Gcbhardt et al. (2000)). This makes it unlikely that our 
results are severely biased by the use of spherical dynam- 
ical models. As discussed in K+2000, the mean intrinsic 
short-to-long axis ratio for our sample of elliptical galaxies 
is< c/a>~ 0.79. Because these are luminous galaxies, ro- 
tation will not contribute substantially to any flattening, 
so except for the flattest galaxies the bias introduced by 
face-on circular orbits will be small. On the other hand, 
possibly embedded face-on disks are likely to be less than 
0.1 — 0.3i? e in size on statistical grounds (Mchlcrt et al. 



(1998)). In three galaxies of our sample, such disks are 
known and extend to 6, 7, and 8 arcsec in NGC 4472, 
NGC 4494, NGC 7626, respectively. In these very inner 
regions, these disks might cause the anisotropy to be over- 
estimated by ~ 0.2, but they will not affect the globally 
averaged results significantly. 

o q ..max _ 
1.6. V c -dO.l 

We have found that both the circular velocity curves 
and the anisotropy profiles of elliptical galaxies are sur- 
prisingly similar. In zeroth order the CVCs can thus be 
characterized by two scaling constants, the effective radius 
and a velocity scale. In Fig. 5 we have used the maximum 
circular velocity; but in fact it should not matter which 
velocity is used to set the scale. In particular, we would 
expect that a suitably defined central velocity dispersion 
could equally be used. A little care is needed, however, 
since the measured central velocity dispersion may be in- 
fluenced by the gravitational field of a central black hole 
and by the resolution of the kinematic data. We therefore 
use an average central velocity dispersion cto.i, defined as 
the square root of the average of all measured of = <r 2 (Ri) 
inside 0.1i? e or 3", whichever is the larger of the two radii. 
This would be the inner luminosity-averaged rms veloc- 
ity dispersion if the surface brightness profile were exactly 
proportional to 

Figure 6 shows a plot of cr .i versus the maximum cir- 
cular velocity for all elliptical galaxies in our sample but 



NGC 4486B. The error-bars for 



are taken from the 



allowed model range as before. For the error-bars in the 
(To.i we have taken the larger of 



Aq-q.i 
o"o.i 



2y/W 



1<1 



(1) 
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Fig. 4. — Velocity anisotropy /3 as a function of radius. Upper panel: for the galaxies with extended data (EK subsample), lower 
panel: for the galaxies from the BSG subsample. 
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Fig. 5. — Average velocity anisotropy /3, calculated over the range 0.1i? e to 17? e or -R max if ^?max < Re, versus maximum circular 
velocity. Solid squares: galaxies with extended data, EK subsample; open squares: galaxies from BSG subsample; see text and 
K+2000. The error-bars plotted for (3 show the standard deviation from the mean in the radial range used in the averaging. 



and a weighted observational error in R < 0.1R e , Aa bs = 

(^Act^ 2 ) . For some galaxies the actual error of cto.i 
may be smaller than this conservative estimate from the 
standard deviation of the of, but we have used this be- 



cause several galaxies have velocity dispersion gradients in 
the central 0.1i? e . Fig. 6 shows a very good correlation; 
the slope is 1.062±0.058 for a reduced \ 2 = 0.87, as deter- 
mined by the routine fitexy of Press et al. ((1992)). This 
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Fig. 6. — Top: the correlation between the averaged central velocity dispersion and the inferred maximum circular velocity. 
Bottom: the residuals from this correlation plotted versus the mean orbital anisotropy defined in Section 2.2. Symbols as in Fig. 5. 



routine fits a straight line to the datapoints by minimiz- 
ing a x 2 function which involves the errors in both the 
x- and y-variables. Thus the fitted relation is completely 
consistent with linear; explicitly, 



(to.1 - 0.66< 



(2) 



The lower panel of Figure 6 shows that the residuals 
from this relation are correlated with the mean orbital 
anisotropy defined in §2.2 (the correlation would be some- 
what stronger had we plotted it in terms of the maxi- 
mum anisotropy). Differences in dynamical structure thus 
cause some scatter in the relation between ao.i and t>™ ax 
but their influence on this and other global correlations 
is small. In this sense the dynamical structure of these 
ellipticals is indeed very uniform. 

3. SCALING RELATIONS 

3.1. Faber-. Jackson relation 

The well-known relation between the total magnitude 
and central velocity dispersion for elliptical galaxies (Faber 
and Jackson (1976)) is shown in the top panel of Fig- 
ure 7 for the K+2000 sample of ellipticals, with oo.i in 
place of the central dispersion. In this plot and all sub- 
sequent similar plots is the twice the total integrated 
B-band absolute magnitude within R e (column 7 of Ta- 
ble 4 in K+2000). The errors in these B l T magnitudes 
are set to 0.3 mag. This accounts for observational er- 
rors and uncertainties in the distances. The latter are 
estimated as approximately 0.25 mag from the intrinsic 
scatter of the Dn-u relation, the depth of clusters (in case 
of group/cluster distances), and a comparison of various 
distance determinations by Tonry et al. ((1997)). The fit- 
ted slope of the Faber-Jackson relation for our ellipticals 
(again using the errors in both variables) is —4.89 ± 1.12, 



i.e., L B cx (cto.i) 1 ' 96 * - 45 - The fit excludes NGC 4486B; 
this galaxy is a close companion of M87 and its low lumi- 
nosity for its high circular velocity suggests that the galaxy 
may be tidally disturbed. The uncertainty in the slope 
is determined with rescaled errors such that the reduced 
X 2 = 1.0. The shallow slope is consistent within our er- 
rors with the slope determined from the data of Faber et al. 
(1989), which give Lb oc cr 2 -6i±o.08 usm g their R-distances. 
For comparison, the K-band slope is significantly steeper 
(Pahre et al. (1998)). 

3.2. Tully- Fisher relation 

The relation between the circular velocity and total lu- 
minosity, known as Tully-Fishcr ((1977), TF) relation, is 
observationally well-established for spiral galaxies. From 
the K+2000 dynamical analysis, we now know circular ve- 
locities also for elliptical galaxies, and can thus investigate 
whether ellipticals also follow a TF relation. We use the 
»™ x variable and errors from §2 and integrated total B l T 
as before, for the 21 galaxies analyzed in K+2000. These 
data points are shown in the middle and lower panels of 
Figure 7. From a least-square fit in both variables we ob- 
tain a B-band TF-slope of —5.92 ± 1.21, corresponding to 

(3) 



Lb oc (v c ) 



max\2.37±0.48 



again excluding NGC 4486B and using rescaled errors to 
estimate the uncertainty. The derived slope is consistent 
with the Faber-Jackson relation and the v™ ax -ao.i relation 
discussed in §2.3. The zero-point of the relation is given 

by 



where Lu = Lb/10 
iJo/eSkms-iMpc" 1 



1 K.65 L &:B 



this implies v" 



and h , 65 
■* ~ 303 kms" 
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Fig. 7. — Top panel: Faber- Jackson relation for our sample of elliptical galaxies. Full squares and open squares show galaxies 
from the EK and BSG subsamples of K+2000. cro.i is the averaged velocity dispersion inside 0.1R e (see text); this is not available 
for NGC 4486B. The full line shows a least-square fit. Middle panel: Tully-Fisher relation for the elliptical galaxies from K+2000 
(same point styles) and SOs from Neistein et al. (1999) (stars). Ellipticals and SOs form a smooth sequence; the slope of the 
elliptical galaxies alone (full line, fit excludes NGC 4486B) is marginally shallower than the slope for the combined sample (dashed 
line). Lower panel: Comparison of elliptical and spiral galaxy Tully-Fisher relations. Ellipticals: Data points and least-square 
fit line repeat the B-band relation. The upper heavy line shows this relation shifted to the R-band, using the colour-magnitude 
relation given in the text. Spirals: Cepheid-calibrated B-band Tully-Fisher relations from Federspiel et al. ((1998), short-dashed) 
and Sakai et al. ((2000), lower long-dashed), Cepheid-calibrated R-band relation from Sakai et al. ((2000), upper long-dashed), 
and r-band relation from Courteau & Rix ((1999), dotted) shifted to the B-band using the colour-magnitude relation for their 
spirals. 
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for an L*-elliptical galaxy, using a corrected Mg(B^) = 
-20.8 (L* B (B°) = 1O 1O - 5 L B , ) for H = GSkms- 1 Mpc" 1 
from Fukugita & Turner ((1991)). The corresponding 
a* A ~ ^kms" 1 . 

We have also included the SOs from Ncistcin et al. 
((1999)) in the TF plot. We took the circular velocities 
corrected for asymmetric drift from their Table 1, column 
13, but used the B-band luminosities from the RC3 (de 
Vaucouleurs et al. (1991)) rescaled to our distance scale 
(using distances from Faber et al. ((1989)) where available: 
NGC 584, 1052, 2768, 3115, 4649, and for the remaining 
SO galaxies redshifts with respect to the CMB frame and 
Ho = 65 km/s/Mpc). The SO's join smoothly with the 
ellipticals in the TF plot. The fitted slope for both sam- 
ples together is slightly steeper (—7.09 ±0.91) than for the 
ellipticals alone, but the two slopes are consistent within 
their errors. 

The lower panel of Fig. 7 shows the comparison of the 
elliptical and spiral galaxy TF relations. The data points 
and the heavy least-square fit line repeat the B-band ellip- 
tical galaxy relation. The two lower dashed lines show the 
Cepheid-calibrated B-band Tully-Fisher relations for spi- 
ral galaxies from Federspiel et al. ((1998)) and Sakai et al. 
((2000)), which are in mutual agreement to a few tenths 
of a magnitude over the range of interest here. The TF re- 
lations of these authors are given in terms of HI linewidth 
measured at 20% peak flux and were converted to u™ ax 
by using Fig. 18 of Rubin et al. ((1999)), which shows 
that the inclination-corrected W20 = 2.0w™ ax to within 
the errors. The dotted line shows the r-band relation from 
Courteau & Rix ((1999)), converted to the B-band using 
the colour magnitude relation they give for their sample 



of spirals, and converted to our Hq = 65 kms -1 Mpc -1 
distance scale. Their velocities V2.2 were used directly as 
v™. The average of the spiral galaxy slopes is somewhat 
steeper than the elliptical galaxy relation, but given the 
large scatter in the latter and between the spiral galaxy 
slopes, this difference is marginal. 

Federspiel et al. ((1998)) applied their calibration to a 
complete sample of Virgo cluster spiral galaxies to derive 
a mean cluster distance modulus of (m-M) = 31.58 ±0.24 
mag. According to Faber et al. ((1989)), the Virgo clus- 
ter center is at a distance of v = 1333 km s -1 . By us- 
ing relative distances from Faber et al. ((1989)) (the R 
values of their Tables 3, 4) and a Hubble constant of 
H = 65kms -1 Mpc -1 to obtain the absolute magni- 
tudes W T for our elliptical galaxies, we have thus im- 
plicitly assumed a Virgo cluster distance modulus of m- 
M = 51og(1333/65) + 25 = 31.56 mag, which is iden- 
tical to the spiral galaxy mean Virgo distance modulus 
of Federspiel et al. ((1998)). Thus the zero points of the 
Cepheid-calibrated TF relations are directly comparable 
to our elliptical galaxy TF. In fact, this comparison need 
not make any assumption about absolute distances; the 
only assumption made is that the centroid of the Virgo 
cluster spiral galaxy sample coincides with the Virgo clus- 
ter ellipticals as given by Faber et al. ((1989)); the relative 
distances between the ellipticals in our sample are then 
also fixed. Therefore, the offset between elliptical and spi- 
ral galaxies in Fig. 7 can be regarded as one in apparent 
magnitude. Its value depends only on having used the 
correct relative distance of the two systems. 

Taking an average over the different spiral galaxy TF re- 
lations in the plot, we conclude that, at a given circular ve- 
locity, elliptical galaxies are about 1 mag fainter in B than 



Dynamical family properties of giant elliptical galaxies 



9 



O 




2.2 



2.4 



2.6 



log 



2.8 



Fig. 9. — Luminous mass in stars versus maximum circular velocity for our sample of elliptical galaxies. The dotted lines show 
the baryonic mass rotation velocity relations for spiral galaxies from McGaugh et al. ((2000), steeper) and Bell & de Jong ((2000), 
shallower) . The heavy parts of these lines cover the range of the fitted spiral galaxy data. 



spiral galaxies. Put differently, at given luminosity ellipti- 
cal galaxies have higher circular velocities than spirals by 
about 0.2 dex. This difference decreases slightly if instead 
of maximum circular velocities we use the halo velocities 
of our models. To see whether this difference is smaller 
in the redder R-band, where elliptical galaxies should be 
relatively brighter, we have plotted in the lower panel of 
Fig. 7 the R-band TF relation expected from the B-band 
fit and the colour-magnitude relation from the EFAR sam- 
ple (Saglia et al. (1997a)), B-R= -0.030(R + 22.5) + 1.50 
(upper full line), to be compared to the R-band Cepheid- 
calibrated relation from Sakai et al. ((2000)), plotted as 
upper long-dashed line. With respect to the Sakai et al. 
relations in R and B, the offset is indeed ^0.4 mag smaller 
in R than in B at logw™ ax '*, but note the different slopes. 

In Figure 8 we show the residuals from the elliptical 
galaxy TF relation versus the mean velocity anisotropy 
/3mean and the effective radius. There is no obvious corre- 
lation with /3 mea n- On the other hand, the plot against R e 
does show a correlation, which reflects the existence of the 
fundamental plane (see next section). 

In summary, elliptical galaxies follow a Tully-Fishcr re- 
lation, with a B-band slope that is marginally shallower 
than the slope for spiral galaxies, and zero point such 
that an L*-elliptical has a circular velocity u™ ax '* ~ 
300 km s -1 . At given circular velocity, elliptical galaxies 
are about 1 mag fainter in B than spiral galaxies. 

3.3. Baryonic Tully-Fisher relation 

Since the stellar mass-to-light ratios of elliptical galax- 
ies vary significantly with luminosity (see §3.5), it may be 
more revealing to plot luminous mass against circular ve- 
locity than luminosity. The mass of the X-ray emitting 
gas in ellipticals is only a few percent even for luminous 
galaxies (Forman, Jones & Tucker (1985), Sarazin (1997)). 
The stellar mass is thus nearly equal to the total baryonic 
mass. Fig. 9 shows stellar mass M* = M/L centra i x Lb 
versus maximum circular velocity for our sample of ellip- 
tical galaxies, where M/L centla \ is the maximum M/L al- 
lowed by the kinematic data. The steeper dotted line in 



the diagram is the best fit spiral galaxy baryonic Tully- 
Fishcr relation relation from McGaugh et al. ((2000)), in- 
cluding stellar and gaseous mass. Their stellar masses 
were determined from red and NIR luminosities, using 
M I 'Us based on population models (constant star forma- 
tion rate, Salpeter IMF, such that M/L K = 0.8 M Q /L Q , K , 
M/L B = IAMq/Lq.b). According to McGaugh et al. 
((2000)) these model M/L are consistent with maximum 
disk fits for the bright galaxies. The shallower dotted line 
is the spiral galaxy line from de Bell & de Jong ((2000)) . 
This is based on luminosities in several passbands and M/L 
ratios determined from evolution models, which use an 
IMF containing fewer low-mass stars than a Salpeter IMF, 
as suggested by recent observations in the Galaxy (e.g., 
Gould, Bahcall & Flynn 1997, Holtzman et al. (1998)) and 
so that maximum disk M/Ls are not exceeded. Both spiral 
galaxy relations apparently predict very similar baryonic 
masses for the circular velocities of interest here. The thin 
parts of both lines are extrapolated from the thick parts 
with spiral galaxy data. 

In Fig. 9 elliptical galaxies fall below the spiral galaxy 
lines by a factor of about two in the mean. Note that the 
offset to the spiral galaxies cannot be explained by lumi- 
nosity errors for our ellipticals alone (Saglia et al. (1997b)). 
The diagram is changed little if instead of t>™ ax we use the 
on average slightly lower circular velocities at lR e for the 
elliptical galaxy points. Note, however, that there is only 
a partial overlap in velocity for the spiral and elliptical 
galaxy samples. With the present data the case for el- 
liptical galaxies having indeed lower baryonic mass than 
spiral galaxies of the same circular velocity is therefore not 
entirely clear. 

It is noteworthy, however, that the baryonic masses of 
elliptical galaxies from dynamics are if anything slightly 
lower than the baryonic masses of spiral galaxies from re- 
alistic stellar population models, in the region where both 
distributions overlap. This suggests that the underlying 
assumption, that ellipticals are described by maximum 
stellar mass models is approximately correct and, hence, 
that elliptical galaxy halos have indeed fairly flat cores: in 
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Fig. 10. — Upper panel: B-band Fundamental Plane for the sample of elliptical galaxies from K+2000. The parametrization of 
the FP is taken from j0rgensen et al. ((1996)). The line shows the best fit with fixed slope 1.0. Lower panel: residuals of the fit 
plotted against mean anisotropy (3. Symbols as in Fig. 5. 



hierarchical models in which ellipticals form through merg- 
ing, a continuity between spirals and ellipticals would be 
expected. Near- maximal M/L in ellipticals are in line with 
results for the Milky Way (Gerhard (1999)) and barred 
galaxies (Debattista & Sellwood (1998), Weiner, Sellwood 
& Williams (2000)), where independent dynamical con- 
straints on the luminous mass favour near-maximal disks, 
while the situation is less clear for luminous unbarred 
spirals (see, e.g.,. Athanassoula, Bosma & Papaioannou 
(1987), Courteau & Rix (1999), Salucci & Persic (1999), 
Bell & de Jong (2000)). 

3.4. Fundamental plane 

In the three-dimensional space defined by central veloc- 
ity dispersion, effective surface brightness (I e ), and effec- 
tive radius (R e ) elliptical galaxies fall on a "fundamental 
plane" (FP; Dressier et al. (1987), Djorgovski & Davis 
(1987)). The existence of the FP is thought to be a conse- 
quence of the virial theorem together with a systematic re- 
lation of M/L on luminosity (Dressier et al. (1987), Faber 
et al. (1987)). Figure 10 shows the FP projection after 
J0rgensen et al. ((1996)) for our sample of ellipticals. Ef- 
fective radii R e are taken from Table 3 of K+2000, effective 
surface brightnesses I e are computed from the (corrected) 
SB-profiles integrated to lR e (Table 4 of K+2000), and 
for the central velocity dispersion we have used the cto.i 
defined in §2.3. The best-fitting line of slope 1.0 is shown 
along with the data in the top panel of Fig. 10. The resid- 
uals with respect to the best fitting line arc plotted in the 
bottom panel against the velocity anisotropy /3 mea n- As 
Fig. 10 shows, these galaxies follow the FP well and the 
residuals are not correlated with /3 me an- The rms scatter in 
log R e is 0.084, excluding NGC 4486B. If we fit the slope 



as well, it is 0.924+ 0.069 and the scatter becomes 0.082. 

In the previous section, we have remarked on the correla- 
tion of the residuals from the elliptical galaxy TF relation 
with R e . This implies that also in the «™ ax — R e — B l T 
space elliptical galaxies fall on a fundamental plane. Be- 
cause «™ ax is very tightly correlated with cto.i (Fig- 6), the 
resulting plane is very similar to the standard FP. 

3.5. Luminous mass-luminosity relation 

From the variables (R e , <7o.i, I e ) one may define a lu- 
minosity L = c^IeRl and mass Mpp — 3cMR e o~o.i/G. 
As first pointed out by Faber et al. ((1987)), the virial 
theorem would then predict I e oc <Tq iRg , provided the 
mass-to-light ratio and the structure constants Cl and Cm 
are identical for all galaxies. The tilt of the FP with re- 
spect to the virial relation therefore implies either a lu- 
minosity dependence of M/L or some deviation of the 
family of elliptical galaxies from homology. In the for- 
mer case and if the FP is R e oc o-§\I^, we would ex- 
pect M/L OC L(2-«)/2a 7( (- 2 -«-W2a_ p or thfi p resen t 

sample a — 1.15 and (3 — —0.758, so that one obtains 

M/L OC i0.37j-0.05_ 

From our dynamical models we have derived central and 
cumulative B-band mass-to-light ratios from observed SB- 
profiles and kinematics. The central M/Lb corresponds 
to the mass-to-light ratio of the stellar population, on the 
assumption of a maximum stellar mass (minimum halo) 
model in which the luminous stars provide as much mass 
in the central parts as is allowed by the kinematic data. 
These M/Lb values have no residual correlation with 
anisotropy (Fig. 11); this is expected as the anisotropy 
was taken into account in the modelling. 
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Fig. 11. — The central B-band M/L versus mean velocity anisotropy (3. No correlation is found. 
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Fig. 12. — Ratio of the cumulative mass M(R e ), luminous and dark, to FP mass 3-R e oo/G, versus circular velocity at R e . 
Symbols as in Fig. 5. 



Figure 12 shows the ratio of the dynamically determined 
cumulative mass at lR e , including any dark mass that 
may be necessary up to this radius, and the FP mass, 
Mpp = 3R e aQ/G, for galaxies with data to beyond R e . 
This ratio measures the structure constant cm- From the 
figure it is clear that there is no systematic trend of cm 
with circular velocity or mass. Because both photometry 
and dynamical anisotropy were taken into account in the 
modelling, this shows that the dynamical non-homology 
mechanism proposed by Graham & Colless ((1997)) is not 
the main cause for the tilt of the FP. On the contrary, 
the dynamical structure of elliptical galaxies is remarkably 
uniform (Fig. 6), and no correlation was found between 
M/L and anisotropy (Fig. 11). 

Figure 13 shows the dynamically determined mass-to- 
light ratios against luminosity. The upper panel of Fig. 13 
shows the central M/Lb, i.e. the inferred mass-to-light 
ratio of the stellar population in our models. The range 
in central M/Lb is about a factor of three. The lower 
panel shows the cumulative mass-to-light ratio at R e . The 
scatter in these diagrams is comparable to eachother, and 
galaxies with particularly low or high M/Lb correspond 
almost one by one. For the EK sample the scatter is simi- 
lar to the scatter in the FP mass-luminosity relation while 



for the BSG sample it is somewhat larger, as might be 
expected. 

Contrary to the range in central Mj Lb, which is fixed by 
the dynamical models, the actual slope of the M/L — L re- 
lation also depends on the extrapolation used to derive the 
total luminosity; the K+2000 values were derived from the 
photometry using a de Vaucouleurs law. The fitted slopes 
are slightly steeper than the predicted FP relation. For 
the EK sample they are for M/L centra i: 0.57 ± 0.11 when 
the errors are rescaled to \ 2 = 1-0, and for M/L(R e ): 
0.59 ±0.09 with \ 2 = 0.7. For the total sample (excluding 
NGC 4486B) they are for M/L centra i: 0.65 ± 0.09 when 
X 2 = 1.0, and for M/L(R e ): 0.67 ± 0.15 when X 2 = 1-0. 
These slopes are different from the predicted FP slope by 
1.5-2.5 times the a of the M/L fit, and they are either 
influenced by small number effects (for the EK sample) or 
by outlying data points (for the full sample). They thus 
appear still consistent with the predicted FP slope. Fig. 6 
and Fig. 12 which demonstrate the dynamical similarity 
of our galaxies also suggest that the difference may not in 
fact be significant. 

From these results we conclude that the M/L variations 
with luminosity indicated by the tilt in the fundamental 
plane are real. The inferred FP M/L ratios for our sam- 
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Fig. 13. — Mass-to-light ratios as a function of luminosity. The upper panel shows the log of the central B-band M/L versus 
total Lb- The lower panel shows the cumulative M/LsiRe) at one effective radius versus Lb, including any dark mass. Symbols 
are as in Fig. 5. Total Lb and errors as in Fig. 7. The errors on the mass-to-light ratios correspond to one-half the range in the 
M/L-profiles of K+2000, as determined by dynamical models bounding the respective confidence interval for each galaxy. In each 
panel two least-square fits are shown, for all galaxies in the sample but NGC 4486B (dashed lines), and for the EK subsample 
only (full lines). The slopes of these two fits are not significantly different. 



pie of ellipticals correspond to the dynamically measured 
M/L values for the stellar populations, assuming a mini- 
mal halo, and they are not due to a non-homologous dy- 
namical structure changing gradually with luminosity. It 
is still possible, however, that photometric non-homology 
influences the slope of the M / L — L-relation. The trend of 
M/L with L is also not caused by an increasing fraction of 
dark matter as luminosity increases, unless two thirds of 
the luminous mass in the cores of the most massive galax- 
ies is dark matter, and luminous elliptical galaxies arc then 
significantly short of baryons compared to spirals (see §3.3 
and Fig. 9). The change of M/L with L is therefore most 
likely due to the population itself. 

4. CENTRAL M/L AND STELLAR POPULATIONS 

The question of what causes the variation with luminos- 
ity of the M/Ls observed in ellipticals has been discussed 
at length in the literature, without any fully satisfactory 
answer. Systematic deviations from homology have been 
argued to play a role (Pahre et al. (1998)), but as discussed 
above this is not a viable solution. It is possible to fit the 
trend observed in the B-band as a metallicity sequence 
of an old stellar population (Maraston (1999)). However, 
the M/Ls in the K band which this model would predict 
are independent of metallicity, and therefore no correla- 
tion with K-band luminosity should be expected, contrary 
to what is observed (Pahre et al. (1998)). Forbes et al. 
((1998)) and Forbes & Ponman ((1999)) analyse the ages 
determined for 88 galaxies and conclude that the observed 
correlation between age and luminosity is far too weak to 



explain the observed M/L trend with luminosity in any 
band. 

Here we explore whether stellar population models can 
reproduce the central B-band M/Lb values derived from 
our (minimum halo) dynamical models. We compare these 
M/Lb values with the predictions of stellar populations 
models in Fig. 14. We use the stellar population models of 
Maraston ((1998)), as employed in Saglia et al. ((2000b)), 
to interpolate the M/Lb expected for a simple stellar pop- 
ulation of given metallicity and age. These models take 
into account stellar evolution mass loss, i.e., they include 
in the mass budget the masses of stellar remnants, but not 
the mass losses of the progenitor stars. This is different 
from what is done by, e.g., Worthey ((1994)) or Bruzual 

6 Chariot ((1996)), where the total initial mass is con- 
served. In addition to the classical Salpeter IMF (with 
power law exponent 7 = —2.35 independent of stellar 
mass), we consider two additional choices for the IMF: the 
recent comprehensive determination of Kroupa ((2000)), 
indicating that the low mass stars could be less numerous 
(his eq. (2): 7 = — 1.3 for to < O.5to , 7 = —2.3 at larger 
masses), and the more extreme IMF of Gould, Bahcall 
and Flynn ((1997), hereafter GBF), where a flatter slope 
at lower masses is suggested (after correction for binaries, 

7 = -0.9 for m < 0.6M Q , 7 = -2.21 for 0.6 < m < 1 and 
7 = -2.35 for to > IM Q ). For all three IMFs we use a 
lower stellar mass cut off of 0.1M Q . 

As input we use ages and metallicities as derived 
by Kobayashi & Arimoto ((1999)), Terlevich & Forbes 
((2000)) and Trager et al. ((2000a)). Such data are avail- 
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Table 1 

The M/L b values of the stellar populations in solar, units, using the models of Maraston ((1998)). Col. 1 

GIVES THE GALAXY NAME, COL. 2 THE M/Lb OF THE LUMINOUS COMPONENT FROM OUR DYNAMICAL ANALYSIS, USING 

Hq = 65kms _1 Mpc -1 , Col. 3 the M/Lb derived using the metallicities of Kobayashi & Arimoto ((1999), KA) 
with 15 Gyr age, Col. 4 using ages and metallicities from Terlevich & Forbes ((2000), TF), and Col. 5 using 

AGES AND METALLICITIES INSIDE (Re/ 8) OF TRAGER ET AL. ((2000a), (2000b), TFWG). COLUMNS 3, 4 AND 5 ARE 

computed for a Salpeter IMF. Cols. 6, 7 and 8 repeat Cols. 3, 4 and 5, but for the Kroupa ((2000)) IMF. Cols. 

9, 10 AND 11 ARE THE SAME FOR THE GBF IMF. 
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able from at least one study for 17 out of 21 of our galaxies. 

Kobayashi & Arimoto ((1999)) derived mean metallici- 
ties inside li? e , considering line index gradients. We use 
values as given in their Table 2. These are derived from 
the Mgi line, assuming an age of 17 Gyr. The M/Lb 
values shown in Fig. 14 are determined by reducing their 
metallicities by 0.2 dex (to correct for the Mg/Fe over- 
abundance) and for an age of 15 Gyr. 

Terlevich & Forbes ((2000)) compiled a catalogue of high 
quality absorption line measurements for galaxies and de- 
rived separate age and metallicity estimates using Worthey 
((1994)) models. 

Trager et al. ((2000a)) determined ages and metallicities 
after applying a correction for the Mg/Fe overabundance 
to the line indices. We use their ages and metallicities for 
the inner R e /8 datapoints, and average over the 4 models 
considered by these authors, for all galaxies except NGC 
1399. For this galaxy we use the new age and metallicity 
determination by Trager et al. ((2000b)); for the other 
galaxies the new determinations agree with the previous 
values within the errors. 

Table 1 lists the resulting stellar population M/Lb val- 
ues with each of the three IMFs. Fig. 14 shows the 
comparison of the stellar M/Ls determined with Kroupa's 
((2000)) IMF with the dynamical M/L B s, rescalcd to 
Hq = 50kms _1 Mpc" 1 which gives the best overall agree- 
ment. The galaxy- by-galaxy comparison after this rescal- 
ing is reasonable within the errors. Only two objects 



are particularly deviant, NGC 315, where the dynamical 
value is a factor 2 larger than the stellar population es- 
timates based on Terlevich & Forbes ((2000)) and Trager 
et al. ((2000a)), but is in agreement with the value ob- 
tained from Kobayashi & Arimoto ((1999)), and NGC 
1399, where M/Lb (dyn) is again a factor 2 larger than the 
stellar M/L B based on Terlevich & Forbes ((2000)), but 
agrees with the value derived from Trager et al. ((2000b)). 
This discrepancy is due to the rather low ages (5 Gyr) 
assigned there, which is also in conflict with the determi- 
nations of Maraston and Thomas ((2000)). 

Similar plots are obtained for the Salpeter and the 
GBF IMFs, when distances are scaled to an optimal 
Hq = 75kms~ 1 Mpc~ 1 for the former and Hq = 
40 km s -1 Mpc" 1 for the latter case. These diagrams 
show that the dynamical M/L B obtained with Hq = 
65kms~ 1 Mpc~ 1 are « 30% larger than the stellar M/Ls 
with the Kroupa ((2000)) IMF, w 60% larger with the 
GBF IMF, and w 20% smaller than the stellar M/Ls 
with the Salpeter IMF. Note that an IMF flatter than 
Salpeter for m > 1M produces stellar M/Ls larger than 
the Salpeter values. 

We conclude that our dynamical Mj L B values, based 
on models maximizing the contribution of the luminous 
component, are compatible with those predicted by stel- 
lar population models, within the uncertainties in the dis- 
tance scale and the poorly known fraction of low-mass 
stars present in giant ellipticals. Only in the case that 
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Fig. 14. — The comparison between the dynamical estimates of the central B-band M/L (transformed to H = SOkms" 1 Mpc -1 ) 
with the predictions of stellar population models of Maraston ((1998)), using Kroupa's ((2000)) IMF. Metallicities and ages used 
in the models are taken from Kobayashi & Arimoto ((1999), filled circles), Terlevich & Forbes ((2000), open circles), Trager et al. 
((2000a), (2000b), filled squares). 



(i) an IMF as fiat at low stellar masses as that of GBF 
is applicable to our elliptical galaxies, and simultaneously 

(ii) a short distance scale (H ~ 80 kms" 1 Mpc -1 ) turns 
out to be correct, would we have underestimated the lu- 
minous masses by as much as a factor w 2 by making the 
assumption of maximum stellar, or minimum dark halo 
mass. For lower values of Hq and/or the other IMFs in- 
vestigated the difference is smaller. For comparison, recall 
that to explain the tilt of the FP as due to an increasing 
fraction of dark matter with luminosity, we would require 
the stellar mass in the most luminous ellipticals to be only 
one third of the inferred dynamical mass. Furthermore, 
the ratio of M/Ls(dyn)/M/Ls (stellar) for this sample of 
elliptical galaxies does not correlate with luminosity. 

Together with the results of Section 3 this suggests that 
the tilt of the FP is a stellar population effect. Most likely 
the main driver is metallicity, but then a secondary pop- 
ulation effect is needed to explain the K-band tilt (Pahre 
et al. (1998)). A larger sample of galaxies with detailed 
dynamical modeling also of the K-band profiles will be 
needed to decide whether this is the final answer to the 
problem of the tilt of the Fundamental Plane. 

5. MASS-TO-LIGHT RATIOS AND MINIMUM HALO 
PROPERTIES 

In this section we first show that the dynamical mod- 
els imply significant radial M/L variations in elliptical 
galaxies. Because the implied local mass-to-light ratios 
are large, we argue that these are the signature of dark 
matter halos. Based on the result of the previous section, 
that the dynamically determined central Mj L b values for 
our ellipticals are consistent with the stellar population 



M/Lbs, we use the dynamical models to delineate the de- 
composition of the elliptical galaxy circular velocity curves 
(CVCs) into luminous and dark contributions, and finally 
study the properties of the implied minimum halos. 

5.1. Global and local mass-to-light ratios 

Figure 15 shows the radial profiles of cumulative 
M(r)/ Lb{t) for all galaxies in the two subsamples, as de- 
rived from the respective "best model" CVC (Fig. 1) and 
luminosity profile. There is considerable variety in these 
cumulative M/is-profiles; even within the galaxies with 
extended kinematics, the ratio of cumulative M/L b at the 
outer data boundary to the central value ranges from con- 
sistent with 1 to < 2 (see also Table 7 of K+2000). This 
suggests a corresponding spread in the efficiency of dis- 
sipational segregation and angular momentum loss of the 
stellar and gaseous component in the dark matter halo 
during the formation process. 

Even though the gradients in cumulative M/Lb are 
modest, those in the local mass-to- light ratio p(r)/ js(r) 
are not. (Here p(r) and js(r) are the inferred mass and lu- 
minosity densities). Because these local p/j-ratios are less 
certain than the cumulative M/L^-values we show these 
in Figure 16 only for the best cases from the EK-subsamplc 
with extended kinematics. As the figure shows, these lo- 
cal M/Lbs become large ^ 20 — 30 in the modelled outer 
parts of these galaxies. This, and corroborating evidence 
from X-ray data (e.g., Matsushita et al. (1998), Loewen- 
stein & White (1999)) argues strongly that the measured 
M/L variations are not due to a slow outward change of 
the stellar population, but instead imply dark matter ha- 
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Fig. 15. — Cumulative B-band M(r)/Ls(r) as function of radius. Upper panel: EK galaxies, lower panel: BSG galaxies. 
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Fig. 16. — Local B-band M/Lb as function of radius for the four most reliable galaxies from the EK sample: NGC 1399 (full 
line), NGC 7626 (dot-dashed), NGC 7507 (long-dashed), and NGC 2434 (short-dashed), over the range of the modelled kinematic 
data. 



los similar to those inferred in spiral galaxies, where the 
component contributing most of the mass, baryonic or not, 
is very different from a normal stellar population. 

5.2. Circular velocity curve decomposition 

Thus, as for spiral galaxies, it is of interest to analyze 
the relative contributions of the luminous and dark mat- 
ter components to the CVCs in Fig. 1. In doing this we 
assume that the M/Lb of the luminous component is con- 
stant with radius and has the maximal value allowed by 
the kinematic data, providing nearly all the mass in the 
center. The dynamical models used by K+2000 to an- 



alyze the kinematic data were built on this assumption, 
but once the CVCs are determined, the luminous compo- 
nent could in principle be assigned less mass a posteriori. 
The discussion in §4 has shown that the maximum central 
M/Lb values determined for our ellipticals by K+2000 are 
consistent with the M/Lb values expected for the stellar 
population of these galaxies, within the uncertainties in 
the distance scale and the lower-mass IMF. Fig. 9 has also 
shown that even with maximal M/Lb elliptical galaxies 
have if anything slightly lower baryonic mass than spiral 
galaxies of the same circular velocity. 

Figure 17 shows the total circular velocity curves from 
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Fig. 17. — Total circular velocity curves (solid lines), and contributions of the luminous matter (dotted line) and the dark halo 
(dashed). Note that the curves are extrapolated to 5R e . 



Fig. 1 for all galaxies in the high-quality EK-sample, and 
the respective contributions from the maximum luminous 
and corresponding minimum dark halo components. At 
li? e , the halo contributes between 1/4 to 2/3 of the cir- 
cular velocity in these "best" dynamical models, corre- 
sponding to between 10 — 40% of the integrated mass. All 
curves are plotted to radii of 5i? e , extrapolating the mod- 
els beyond the radial range of the data. The case of NGC 
315 (where the outer rise of the CVC is due to modelling 
problems, see K+2000) shows that such extrapolation can 
lead to large errors. However, in cases where X-ray data 
or planetary or globular cluster velocities were available, 
the "best" models of K+2000 matched the independent 
mass estimates from these outer data very well. Figure 17 
shows that within the framework of the models, luminous 
and dark matter reach equal interior mass at <~ 2 — 4i? e , 
and at 5R e the halo is predicted to dominate in all models 
except in one case. As in spiral galaxies, the combined ro- 
tation curve is flatter than that for the individual compo- 
nents ("conspiracy"); this is already seen within the radial 
range of the kinematic data. 

The last panel of Figure 17 repeats the CVC decompo- 
sition for NGC 2434 with a luminous M/Lb = 3.4 instead 
of the maximum M/Lb = 6.0. This is the average of 
the population values obtained for this galaxy with the 
Kroupa and GBF IMF's (see Table 1). The lower M/L B 
leads to a significantly denser halo; the decomposition in 
this case is comparable to those of Rix et al. (1997)) with 
a Navarro, Frcnk & White (NFW, (1996)) halo mass dis- 
tribution. This example suggests (i) that the kinematic 
data for the K+2000 sample could presumably have been 
fit also by using these halo models, and (ii) that the re- 
sulting lower M/Lb values for the luminous components 



would be near the lower end of the range consistent with 
the stellar population models of Section 4. However, in 
this case, ellipticals would move further down in Fig. 9: 
in the models of Rix et al. ((1997)) with a NFW halo for 
NGC 2434 the M/L for the luminous component is an- 
other (for the same distance) factor of 1.5 smaller than 
the maximum M/L found by K+2000, which had already 
placed NGC 2434 by a factor of 4 below the Bell & de 
Jong ((2000)) and McGaugh et al. ((2000)) lines in Fig. 9 
(at logu™ ax = 2.55). 

In a recent study, Loewenstein & White ((1999)) con- 
sidered the implications of the observed X-ray tempera- 
ture - velocity dispersion relation for the mass distribu- 
tion and dark halos in elliptical galaxies. They found that 
the total (luminous plus dark) M/Ly at 6R e is nearly in- 
dependent of galaxy luminosity, with value M/Lv(QR e ) — 
25ft-8o Mq/Lqy, so that the ratio of dark to luminous mat- 
ter decreases with luminosity. Their M/L(6R e ) converted 
to the B-band and rescaled to is consistent with the 
Mj L B (6R e ) predicted by our models for L B isj L* galax- 
ies, while for our fainter galaxies the dynamical models 
predict ~ 0.2 dex lower mass per luminosity than given 
by Loewenstein & White ((1999)), in all cases, however, 
extrapolating the models beyond the radial range of the 
kinematic data. The dark matter mass fraction within R e 
found from the X-ray analysis (> 20%) is in good agree- 
ment with our result quoted above. 

5.3. Dark halo parameters for ellipticals 

In this section we investigate scaling laws for the halo 
densities and halo core radii of elliptical galaxies based on 
the "minimum halo models" of K+2000. We also compare 
the inferred halo properties with the dark matter halos in 
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Fig. 18. — From top to bottom: (i) Effective radius in kpc versus total luminosity in solar units, (ii) Halo core radius in kpc, (iii) 
halo velocity in kms -1 , (iv) central halo density in Mq/pc 3 , and (v) central halo phase-space density in Mq/ pc 3 /( kms -1 ) 3 , 
all referring to the "best" models selected by K+2000 from the middle of the respective confidence range for all sample galaxies. 
Symbols are as in Fig. 5. The full lines show the least-square fits. Dotted lines show the corresponding relations for spiral galaxies 
given in the text. 



spiral galaxies. 

Figure 18 shows, as a function of galaxy luminosity, the 
effective radius R e , the halo core radius r c the halo ve- 
locity Vh, core density ph, and central phase-space density 
f h , for the respective "best" models of K+2000. All these 
correlate with luminosity, but the tightest correlation is 
that between the effective radius and the luminosity (one 
of the projections of the FP). 

We note that for many of the galaxies in Fig. 18 a con- 
stant M/ Lb model with only luminous mass is within the 
95% confidence range of the K+2000 models. For these 
galaxies the halo radii, circular velocities, and densities 
have therefore large uncertainties in a logarithmic plot 
like Figure 18. However, in cases where X-ray data or 
planetary or globular cluster velocities were available, the 
"best" models of K+2000 matched the independent mass 
estimates from these outer data very well, while the con- 
stant M/L models were usually inconsistent with these 
data. For seven of the sample galaxies, a constant M/L 
model was found inconsistent with the kinematic data 
(K+2000). For these galaxies we have estimated 95% con- 



fidence bounds on the halo parameters from y 2 contour 
plots and plotted them as error bars in Fig. 18. In a few 
of these cases only the halo density is well-determined. 
Within the considerable uncertainties, the 'error bars' are 
consistent with the scatter of the points. Note that some 
of the galaxies with the best evidence for dark halos (NGC 
2434, 7507, 7626) are among the smallest halo core radii 
and largest halo density points with respect to the mean 
fit lines in Fig. 18. By using the mean scaling relations 
from all the "best" halo models of K+2000 for the subse- 
quent discussion, we have therefore not biassed the nor- 
malisations of the best-fit lines, while the best-fit slopes 
are much better determined. 

Figure 18 shows that, in the mean, more luminous galax- 
ies have larger halo core radii, with a slope similar to the 
log(i? e )-log(L)-relation. They also have larger halo circu- 
lar velocities, and lower central densities and phase-space 
densities. As predicted by hierarchical models less massive 
objects are denser. The least-square fits shown in Figure 
18 correspond to the following scaling laws: 

ii e = 11.8L? 1 81 / l ^ 5 kpc (5) 
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r Clh = 13.8LSi 92 ^J 5 kpc (6) 
v h = 397L° 1 fkms- 1 (7) 
p fc = O.O46L7/ 1 1 - O6 ^. 66 M /pc 3 (8) 

fh = 2.0 x lO-^^-'^LsMo/pcVCkms- 1 ) 3 (9) 
where L n = L B /10 n hQ% 5 L QtB , 

* " *4 <10) 

for the employed halo models (see eqs. (2-4) in Gerhard 
et al. (1998)), and the central phase-space density is de- 
fined by 

f h = 2*/*p h /vl (11) 

Eqs. (5) and (6) result in a ratio of R e and halo core radius 
that is approximately constant, 

r c . h /R e = 1.2L?i n . (12) 

Also, from the results of §3.5 and eqs. (5) and (10) one sees 
that the ratio of luminous mass density [M I Lb)L b I R\ 
and dark halo density has little luminosity dependence in 
the mean, although the scatter is large; the density ratio 
varies between 10~ (0 - 5 ~ 2) . 

The dotted lines in Fig. 18 show the scaling relations 
for spiral galaxies from Persic, Salucci & Stel ((1996a), er- 
ratum (1996b); PSS), which are based on the same (min- 
imum) halo models. The relation for spiral galaxy halo 
central densities as given by PSS and rescaled to the dis- 
tance scale used here becomes 

p^O.OO19i n °- 7 4 65 M /pc 3 . (13) 

For luminosities around Lb — 10 11 Lb ;Q ~ 3i», ellipti- 
cal galaxy halos are therefore about 25 times denser than 
spiral galaxy halos of the same Lb, assuming maximum 
stellar mass in both cases. This result agrees well with the 
work of Bertola et al. ((1993)), who used extended HI disks 
around a few elliptical galaxies to constrain their halo mass 
distributions. Most of the factor 25 can be traced back to 
the fact that the CVCs of both spirals and ellipticals are 
approximately flat, and that for an elliptical galaxy pro- 
file, the maximum circular velocity occurs at significantly 
smaller radius in units of R e than for an exponential disk. 
There is an additional factor ~ 2 because of the larger 
M/Lb of ellipticals at given L B - 

According to PSS, the halo core radii of spiral galaxies 
also scale with luminosity when expressed in units of the 
optical radius. We can compare their relation to the case 
of elliptical galaxies as follows: Fitting the parameters R e 
and L for <~ 200 spiral galaxies from the RC3, as given by 
Burstcin et al. ((1997)), results in the least-square fit line 
(rescaled to H = 65kms _1 Mpc~ ) 

flf = 9.0L° 1 53 / l65 1 kpc. (14) 

This is somewhat shallower than the corresponding rela- 
tion for elliptical galaxies, eq. (5), but the R e values of 
spirals and ellipticals are very similar around L — L*. 
Next, we use i?f pt = 3.2Rd = 1.9Rf^ 2 , where i£f pt , Rd 

and Rf , 2 are the optical, scale and half-mass radius for 



an exponential disk, and assume a mean i?f = \.2Ro for 
the distribution between face-on and edge-on. From the 
relation given by PSS we then obtain 

r s cM /R s e = 5.0L° n 2 (15) 

and hence 

r c s ,=45i; i 73 ^ 5 k P c. (16) 

Comparing these relations with eqs. (12) and (6), one sees 
that the minimum halos of spiral galaxies of the same Lb 
and R e have about 4 times larger halo core radii than are 
inferred for elliptical galaxies from our dynamical analysis, 
and that the ratio is only slightly smaller at given Lb if the 
mean Rf [eq. (16)] is used. In view of the uncertainties in 
the transformations, the slopes in these relations appear 
consistent with eachother. 

Because of the luminosity offset in the TF relation, it 
may be more appropriate to compare elliptical and spiral 
galaxy halos at the same baryonic mass or circular veloc- 
ity than at the same luminosity. Because this would mean 
comparing an elliptical galaxy with a spiral of higher lumi- 
nosity, and because of the luminosity dependences in the 
spiral galaxy relations eqs. (16) and (13), this will increase 
the differences found above. Using -y™ ax for the compari- 
son as in Fig. 7, the density ratio Ph/pf thus increases by 
a factor of <~ 2. This may be an overestimate, however, 
as the inferred asymptotic halo velocities Vh are formally 
lower than -y™ ax [by ~ 0.1 dex, compare eqs. (4) and (7)]. 
We note here that the Vh value for L* B predicted by eq. (7) , 
253 km s" 1 , is in agreement with the galaxy-galaxy lensing 
result of Wilson et al. ((2000)) at radii ~ 100 kpc, but em- 
phasize that in our models vh is much more uncertain than 
v max Making the comparison at constant baryonic mass 
would increase the density ratio by a factor less than 1.6, 
from Fig. 9 and eq. 13. Thus we conclude conservatively 
that the halos of elliptical galaxies are at least 25 times 
denser than the halos of spiral galaxies of similar baryonic 
mass or circular velocity. 

These results also suggest that the phase-space densi- 
ties for the minimum halos of elliptical galaxies are higher 
than for those of spiral galaxies of similar baryonic mass. 
From the TF relation spiral galaxies at given Lb have 
circular velocities about 0.2 dex lower than the y™ ax of el- 
liptical galaxies, but perhaps only 0.1 dex lower than the 
Vh values used in eqs. (7) and (9). This is not sufficient 
to compensate their higher densities. Higher phase-space 
densities would rule out stronger adiabatic contraction as 
the explanation for the denser halos in ellipticals. Recall 
also the similar R e of ellipticals and spirals of the same 
Lb- The argument could be circumvented if spiral galax- 
ies had sub-maximal disks and cuspy halos and their halo 
densities had been underestimated significantly by the PSS 
models. However, note that both in the Milky Way (Ger- 
hard (1999)) and in barred galaxies (Dcbattista & Sell- 
wood (1998), Weiner et al. (2000)), where independent 
dynamical constraints on the luminous mass are available, 
the galactic disks are near-maximal, and that the recent 
stellar population models of Bell & de Jong ((2000)) sug- 
gest that high-surface brightness spirals are generally close 
to maximal disks. 

One possible explanation for the much larger densities 
and probably phase-space densities of elliptical galaxy ha- 
los might be that some of the dark matter inferred in the 
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Upper panel: B-band M/L B , normalized by the central M/L B , versus the gravitational acceleration v 2 /r for all 
galaxies from the EK subsample. The "bending upwards" occurs at accelerations a ~ 10~ 9 msec~ 2 , a factor of 10 higher than 
typical in spiral galaxies. Lower panel: Normalized M/Lb versus mean interior density. 



inner regions of elliptical galaxies is baryonic. The evi- 
dence for large M/L-ratios in ellipticals at large radii, of 
order one hundred, from both X-ray (e.g., Matsushita et al. 
(1998)) and weak lensing data (Griffiths et al. (1998), Wil- 
son et al. (2000)), together with the lack of microlensing 
towards the LMC in the Milky Way (Alcock et al. (2000)), 
does however not make this an attractive explanation for 
galactic halos in their entirety; a separate inner baryonic 
dark matter component would be needed. 

An alternative possibility is that most elliptical galaxies 
formed at high rcdshift from progenitors with higher den- 
sities than seen in present-day spiral galaxies. If halo core 
densities are proportional to virial densities, which in turn 
depend on the density of the Universe at the time of col- 
lapse, then the result above implies that elliptical galaxy 
halos have collapsed at redshifts Ze ^ 25 1 / 3 (1 + zs) — 1, 
i.e., ze ^ 5 if the halos of spiral galaxies of similar lu- 
minosity formed at redshifts z$ > 1. Thus our result 
may indicate that giant elliptical galaxies are old, con- 
sistent with evidence from the fundamental plane (van 
Dokkum & Franx (1996), van Dokkum et al. (1998), Ben- 
der et al. (1998)) and line-strength indices (Bender, Ziegler 
& Bruzual (1996)). Unfortunately, while this argument is 
plausible, it is not conclusive until the relation between 
the shallow central halo profiles inferred from observations 
and the cuspy halos predicted from hierarchical collapse 
of dark matter (Navarro, Frenk & White (1996)) is under- 
stood. Moreover, it would appear to be at odds with the 
observation of substantial merging in a moderate z=0.83 
cluster (van Dokkum et al. (1999)), unless the progeni- 
tors had unusually high halo densities also, and also with 
the expectation that some elliptical galaxies should have 



formed recently from mergers of normal spiral galaxies 
(e.g., Schweizer (1998)). 

The plot of halo core radius against luminosity displays 
considerably larger scatter for both subsamples than the 
plot of R e versus L and consequently also the inferred halo 
density shows considerable scatter at a given luminosity. 
What is the origin of this increased scatter? First, it is pos- 
sible that modelling uncertainties contribute to the larger 
scatter in the derived r Ci h- However, we do not think that 
this can be the whole explanation. The three galaxies at 
the upper boundary of the points in the r Ci /j-L-plot which 
have the best kinematic data for their kind, NGC 1399, 
NGC 3379 and NGC 6703, all have log(r ch /R e ) = (0.5 : 
0.6) and show evidence for small if any amount of dark 
matter within the modelled range. On the other hand, 
the best-determined galaxies near the lower boundary with 
the best evidence for additional dark matter, NGC 2434, 
7507, 7626, have log{r ch /R e ) = (-0.15 : -0.35). Fig. 17 
shows that the two groups have significantly different CVC 
shapes for the visible component only: for the first group 
the visible CVC is almost coincident with the dynamically 
inferred total CVC to lR e , whereas for the second group 
there are significant mass discrepancies already at lR e . 
Thus elliptical galaxies at fixed L b appear to have a range 
of luminous matter CVCs and hence dark matter CVCs, 
even though the total CV rotation curves are fairly simi- 
lar. Also, while the most rapidly rotating galaxies in the 
sample, for which we would expect the largest systematic 
errors in the modelling, have predominantly positive resid- 
uals with respect to the least-square line, the non-rotating 
galaxies populate the entire distribution of residuals in- 
cluding the extremes. Thus we believe that most of the 
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scatter in the inferred log(r c h/R e ) is not due to modelling 
effects, but reflects physical differences between the sample 
galaxies. 

Could the amount of dark mass in the centers of prefer- 
entially those galaxies with apparently large r c h have been 
underestimated with our minimum halo models? Then 
we would expect that the derived central M/Lb values 
of these galaxies should be systematically high for their 
luminosity, i.e., we would expect a correlation of positive 
residuals from the M/Lb — L relation in Fig. 13 with posi- 
tive residuals in r c h — L. However, the M/Lb residuals for 
galaxies with large and small halo core radii for their lumi- 
nosities do not show a systematic difference in the present 
sample. 

Thus, we believe the most likely explanation for the 
larger scatter in r c h is that galaxies of similar luminosity 
have different dark matter core radii and central densitites 
determined by the particulars of the merging process in 
which they were made. In this case the scatter in inferred 
halo density (an order of magnitude or more) should per- 
haps reflect mainly the halo densities of the progenitors at 
the time of formation, with the highest (lowest) densities 
corresponding to the earliest (latest) mergers. We have 
attempted to test this by plotting the minimum halo den- 
sities from Fig. 18 versus the population ages of Section 4, 
but no convincing correlation is seen in the present data. 

Figure 19 finally shows mass-to-light ratios M/Lb nor- 
malized by the central value, as a function of accelera- 
tion a = v 2 (r)/r and mean interior density p(< r) = 
3Vc(r) / 4:irGr 2 , for all galaxies in the EK subsample with 
extended data. The estimated uncertainty in a from that 
in v c and distance is of order 50%. Figure 19 shows that 
the "bending upwards" which indicates the onset of dark 
matter takes place at systematically higher accelerations a 
than in spiral galaxies (McGaugh (1999)), by about one or- 
der of magnitude, which is a consequence of the higher halo 
mass densities in elliptical galaxies. This suggests that the 
acceleration scale found in spiral galaxies is not universal 
as required by the modified gravity theory MOND (Mil- 
grom (1983)), i.e., that in elliptical galaxies an additional 
gradient in M/L would be required besides MOND. 

6. CONCLUSIONS 

Based on a uniform dynamical analysis of photomet- 
ric and line-profile shape data for 21 mostly luminous, 
slowly rotating, and nearly round (17 E0/E1, 4 E2) el- 
liptical galaxies by Kronawitter et al. ((2000)), we have 
investigated the dynamical family relations and dark halo 
properties of ellipticals. Our main results are as follows: 

(1) The circular velocity curves (CVCs) of elliptical 
galaxies are flat to within ~ 10% for R > 0.2i? e to at least 
R > 2R e , independent of luminosity. This argues against 
strong luminosity segregation in the dark halo potential. 

(2) Most ellipticals are moderately radially anisotropic, 
with average (3 ~ 0—0.35, again independent of luminosity. 

(3) The dynamical structure of ellipticals is surprisingly 
uniform. The maximum circular velocity is accurately pre- 
dicted by a suitably defined central velocity dispersion. 

(4) Elliptical galaxies follow a Tully-Fishcr (TF) rela- 
tion with marginally shallower slope than spiral galaxies. 
At given circular velocity, they are about 1 mag fainter in 
B and about 0.6 mag in R, and appear to have slightly 
lower baryonic mass than spirals, even for the maximum 



M/Lb allowed by the kinematics. 

(5) The residuals from the TF and Fundamental Plane 
(FP) relations do not correlate with dynamical anisotropy 
p. 

(6) The luminosity dependence of M/L indicated by the 
tilt of the FP corresponds to a real dependence of dy- 
namical M/L on L. The tilt of the FP is therefore not 
due to deviations from homology or a variation of dynam- 
ical anisotropy with L, although the slope of M/L versus 
L could still be influenced by photometric non-homology. 
The tilt can also not be due to an increasing dark matter 
fraction with L, unless (i) the most luminous ellipticals 
have a factor > 3 less baryonic mass than spiral galaxies 
of the same circular velocity, (ii) the range of IMF is larger 
than currently discussed, and (iii) the IMF or some other 
population parameter varies systematically along the lu- 
minosity sequence such as to undo the increase of M/L 
expected from simple stellar population models for more 
metal-rich luminous galaxies. This seems highly unlikely. 

(7) The tilt of the FP is therefore best explained as a 
stellar population effect. Population models show that the 
values and the change with Lb of the maximal dynamical 
M/Lbs are consistent with the stellar population M/LbS 
based on published metallicities and ages within the un- 
certainties of IMF and distance scale. The main driver is 
therefore probably metallicity, and a secondary population 
effect is needed to explain the K-band tilt. 

(8) The population models show that we would have 
underestimated the luminous masses by as much as a fac- 
tor w 2 only if (i) the flattest IMFs at low stellar masses 
discussed for the Milky Way are applicable to our ellipti- 
cal galaxies, and simultaneously (ii) a short distance scale 
(Hq ~ 80 km s _1 Mpc -1 ) turns out to be correct. For lower 
values of H and/or the other IMFs investigated in Section 
4 the difference is smaller. Together with (4) this makes it 
likely that elliptical galaxies have indeed nearly maximal 
M/Lb ratios (minimal halos). 

(9) Despite the uniformly flat CVCs, there is a spread in 
the ratio of the CVCs from luminous and dark matter, i.e., 
in the radial variations of cumulative mass-to-light ratio. 
The sample includes galaxies with no indication for dark 
matter within 2R e , and others where the best dynamical 
models result in local M/Lbs of 20-30 at 2R e . As in spi- 
ral galaxies, the combined rotation curve of the luminous 
and dark matter is flatter than those for the individual 
components ("conspiracy"). 

(10) In models with maximum stellar mass, the dark 
matter contributes ~ 10 — 40% of the mass within R e . 
Our flat rotation curve models, when extrapolated beyond 
the range of kinematic data, predict equal interior mass of 
dark and luminous matter at <~ 2 — 4i? e , consistent with 
results from the X-ray temperature - velocity relation. 

(11) Even in these maximum stellar mass models, the 
halo core densities and phase-space densities are at least 
<~ 25 times larger and the halo core radii <~ 4 times smaller 
than in spiral galaxies of the same circular velocity. Cor- 
respondingly, the increase in M/L sets in at ~ 10 times 
larger acceleration than in spirals. This could imply that 
elliptical galaxy halos collapsed at redshifts z > 5 or that 
some of the dark matter in ellipticals might be baryonic. 
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